Two-species branching annihilating random walks with one offspring 
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We study the effects of hard core (HC) interactions between different species of particles on 
two-species branching annihilating random walks with one offspring(BAW2(l)). The single-species 
model belongs to the directed percolation (DP) universality class. In the BAW2(1) model, a particle 
creates one particle of the same species in its neighborhood with the probability <j(1 —p) and of the 
different species with (1 — cr)(l — p), where p is the hopping probability. Without HC interactions, 
this model always exhibits the DP-type absorbing transition for all a. Even with HC interactions, 
the nature of the phase transitions does not change except near a — 0, where the HC interaction 
destabilizes and completely wipes away the absorbing phase. The model is always active except at 
the annihilation fixed point of zero branching rate (p = 1). Critical behavior near the annihilation 
fixed point is characterized by exponents /3 = v± = 1/2 and vu = 1. 
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During last decades, there have been considerable ef- 
forts to understand nonequilibrium absorbing phase tran- 
sitions from an active phase into an absorbing phase 
which consists of absorbing states M. Once the system 
is trapped into an absorbing state, it cannot escape from 
the state forever. Various one dimensional lattice models 
exhibiting absorbing transitions have been studied and 
most of them turn out to belong to one of two universal- 
ity classes, directed percolation (DP) and directed Ising 
(DI) universality class. While models with the Ising sym- 
metry between absorbing states belong to the DI class, 
models in the DP class have no symmetry between ab- 
sorbing states . 

It is generally accepted that the dimensionality and the 
symmetry between absorbing states play an important 
role in determining the universality classes as in equilib- 
rium critical phenomena. To find out new universality 
classes, it is natural to study models with higher sym- 
metries than the Ising symmetry. To achieve a higher 
symmetry such as the Potts symmetry, one may increase 
the number of absorbing states. However models with 
higher symmetries investigated so far turn out to be al- 
ways active and only critical at the annihilation fixed 
point of zero branching rate S. Recent field theoreti- 
cal study may provide a possible explanation for this j| . 
Although there is no stable absorbing phase, critical be- 
havior near the annihilation fixed point are non-trivial 
and form new universality classes [Ho . 

In systems with more than two equivalent absorbing 
states, there are various kinds of domain walls (or kinks) 
that cannot cross over each other to retain the order 
of absorbing domains. Multiple occupations of domain 
walls at a site is also forbidden. These restrictions may 
be naturally realized by introducing hard core (HC) in- 
teractions between different domain walls. These domain 
wall dynamics can be considered as a multi-species par- 
ticle dynamics with HC interactions. 

Recently, Cardy and Tauber introduced an interest- 
ing multi-species model, A-species branching annihilat- 



ing random walks with two offsprings(Af-BAW(2)), which 
can be solved exactly for all A~ > 1, using renormaliza- 
tion group techniques in bosonic type formulation which 
ignores HC interactions j|. The A r -BAW(2) model is 
a classical stochastic system consisting of N species of 
particles, Ai (i — 1, . . . , AT). Each particle diffuses on a 
d-dimensional lattice with two competing dynamic pro- 
cesses: pair annihilation and branching. Pair annihila- 
tion is allowed only between identical particles (Ai+Ai — > 
0). In the branching process, a particle Ai creates two 
identical particles in its neighborhood (Ai — » Ai + 2Aj), 
with rate a for i = j and rate cr'/(N — 1) for i ^ j. For 
N = 1, this model exhibits a Dl-type absorbing transi- 
tion at finite branching rate . For AT > 1 , this model 
is always active except at the annihilation fixed point of 
zero branching rate. 

Our recent study on the one-dimensional AT-BAW(2) 
model shows that the HC interaction between different 
species of particles drastically changes the universality 
class in a non-trivial way p]. It contradicts the con- 
ventional belief that the HC interaction is irrelevant to 
absorbing-type critical phenomena, because the particle 
density is so low near an absorbing transition that the 
probability of multiple occupations at a site should be too 
small to be significant. However, in multi-species models 
of A^ > 1, the asymptotic density decay near the annihila- 
tion fixed point becomes nontrivial due to the convective 
displacement generated by the HC interaction. 

In this paper, we study two-species branching an- 
nihilating random walks (BAW) with one offspring 
(BAW2(1)) in one dimension. The single-species BAW 
model belongs to the DP universality class • So this 
model is a natural multi-species extension of the DP-type 
systems. We numerically study the effect of the HC in- 
teraction on the absorbing phase transition and find that 
the HC interaction changes the phase diagram consider- 
ably. 

The BAW2(1) model is a stochastic system consisting 
of two species of particles, A and B. Each particle hopes 
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to a near- neighboring site with probability p or creates a 
particle on a near-neighboring site with probability 1 — p. 
The created particle can be either of the same species as 
its parent with probability a or of the different species 
with probability a' = 1 — a. With HC interactions, any 
dynamics resulting in the multiple occupations at a site 
by different species of particles is forbidden. 

We perform dynamic Monte Carlo simulations. With 
the initial condition of two near-neighboring identical 
particles, we measure the survival probability P(t)(the 
probability that the system is still active at time i), the 
number of particles N{t) averaged over all runs, and the 
mean distance of spreading R{t) averaged over the sur- 
viving runs. 

At criticality, these quantities scale algebraically in the 
long time limit as P(t) ~ t~ 5 , N(t) ~ and R(t) ~ i 1/z 
p2[ . The effective exponents defined as 

-8(t)=\og[P(t)/P(t/b)}/logb (1) 

and similarly for r)(t) , 1/ z(t) show the straight line at crit- 
icality and the upward (downward) curvature in the ac- 
tive (absorbing) phase. With b = 5, we estimates asymp- 
totic values of dynamic critical exponents 8, rj, z and 
critical hopping probability p c for various values of a. 

Fig. f shows a — p phase diagram of the BAW2(1) 
model. Without HC interactions, this model always ex- 
hibits the DP-type absorbing transition for all a. For 
example, at a = 0, we estimate the critical probabil- 
ity Pc = 0.785(5) and dynamic scaling exponents 8 — 
0.155(5), rj = 0.32(1) and l/z = 0.640(5) which agree 
well with the DP values 0. 

Even with HC interactions, the nature of the phase 
transitions does not change except near a = 0. How- 
ever, near a = 0, the HC interaction destabilizes and 
completely wipes away the absorbing phase. The model 
becomes always active except at the annihilation fixed 
point of zero branching rate (p c = 1). Numerically, we 
could not pinpoint when the absorbing phase is com- 
pletely squeezed out, due to statistical errors. With 
present simulations, we estimate the upper bound of a as 
<j c ~ 0.05. To identify the scaling behavior near the an- 
nihilation fixed point of the system with the HC interac- 
tion, we analyze the finite-size effects on the steady-state 
particle density p s . 

The scaling behavior near criticality is characterized 
by exponents /3, v±_ and i/u defined as 

P (t) ~t- a , p s ~A^, (2) 

where A = p c — p, £ is the correlation length, r the char- 
acteristic time, p(t) the particle density at time t, and 
p s the steady-state particle density. At the annihilation 
fixed point, the exponents a and z should follow from the 
simple random walk exponents a — l/z and z = 2 Jl3| . 

Using the finite-size scaling theory on the steady-state 
particle density p s |Q 



p s (A,L)=L-^F(ALV^), (3) 

the value of v± is determined by collapsing data of p s 
with P/f± = 1. We measure p s in the steady state, av- 
eraged over 5 x 10 3 ~ 5 x 10 4 independent samples for 
several values of A (5 x 10~ 4 ~ 0.05) and lattice size L 
(2 5 ~ 2 9 ). We find u± = 0.49(4) (Fig. 2). With a = 1/2 
and z — 2, we estimate 

= 1/2, v ± = l/2, i/|| = l. (4) 

The single-species BAW(l) model is one of the sim- 
plest model belonging to the DP class. BAW(l) dynam- 
ics involves the spontaneous annihilation of a particle 
(A — » 0) which is a common feature of DP models. A 
single particle is spontaneously annihilated by the com- 
bination of a branching (creation) and hopping process 
such as A — > A + A — > 0. 

As a decreases, the pair annihilation by diffusion is 
hard to occur, because chance to meet the same species 
of particles decreases. It makes the system more active so 
the critical hopping probability (p c ) must be an increas- 
ing function of a' = 1 — a . This is consistent with our 
numerical phase diagram (Fig. 1). 

The effect of the HC interaction on the phase diagram 
turns out to be rather tricky. For large a region, the 
system tends to form large A- or B-dominated domains. 
The HC interaction accelerates the pair annihilation pro- 
cess because the domain boundary will induce the diffu- 
sion bias to the center of each domain. Therefore the 
system becomes less active and the absorbing phase ex- 
pands with the HC interaction. 

However, for small a region, the above story is com- 
pletely reversed. The system tends to have a high density 
of locally ordered AB pairs. Without the HC interaction, 
most of the pair annihilation processes should occur by 
diffusions across the domain boundary. So, in this case, 
the HC interaction decreases the chance of the pair anni- 
hilation processes and the system becomes more active. 
This implies that the critical lines with and without HC 
interactions should cross each other like in Fig. 1. 

At a = 0, a single particle cannot be spontaneously 
annihilated by the combination of a single branching pro- 
cess and diffusions. It needs at least three branching pro- 
cesses, e.g., A — > AB — » ABA — » ABAB, which can be 
annihilated by successive diffusions in systems without 
HC interactions. So the absorbing phase is still stable 
and the DP nature is maintained. 

However, with HC interactions, the ordered AB pairs 
can not be annihilated by diffusions. A single particle 
has a finite probability to survive asymptotically and the 
entire absorbing phase becomes unstable. The system 
becomes always active except at the annihilation fixed 
point of zero branching rate. 

Critical behavior near the annihilation fixed point can 
be extracted by balancing the particle density changes 
due to branching and pair annihilation processes. Let us 
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consider a particle A created by a B at time t. Annihi- 
lation of the created particle A occurs by encountering 
an independent A particle through diffusions. The time 
scale of this pair annihilation process is given by ordinary 
diffusions, i.e., ~ d 2 where the mean distance between 
particles d is order of pit) -1 . Branching process of the 
parent particle B increases the particle density with the 
time scale Tb ~ (1 — p)^ 1 ■ By balancing these two time 
scales, we expect that the steady-state particle density 
scales as p s ~ (1 — p)^ with (3 = 1/2. This is consistent 
with our numerical finding in Eq.(0). 

The above argument should be valid for general N- 
species BAW(l) models for N > 2 near the annihilation 
fixed point. Preliminary numerical results for N = 3 con- 
firm our predictions. The universality class characterized 
by exponents (3 = v± = 1/2 and u\\ = 1 also includes the 
TV-species BAW(2) model for N > 2 with static branch- 
ing processes, where two offsprings are divided by their 
parent such as 0£?0 — > ABA Q] . It is clear that our above 
argument also applies to this model. 

In summary, we studied the two-species branching an- 
nihilating random walks with one offspring. This model 
exhibits a DP-type absorbing phase transition as in the 
single-species model. With hard core interactions be- 
tween different species of particles, the phase diagram 
changes considerably. Especially, the system becomes al- 
ways active with an annihilation fixed point for very small 
branching rate of the same species. We find that the crit- 
ical behavior near the annihilation fixed point is identical 
to that in the A-BAW(2) model with hard core interac- 
tions and static branching. The A-species generalization 
of the BAW(l) model is currently under study. 
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FIG. 1. The a — p phase diagram for the BAW2(1) model. 
Open and filled circles correspond to the critical points with- 
out and with the hard core interaction respectively. Lines 
between data points are guides to the eyes. 



FIG. 2. Data collapse of p s L 0/l,J - against AL 1/l,± with 
f3/v±_ = 1 for various system size L = 2 5 , . . . , 2 9 for BAW2(1) 
with the hard core interaction at a = 0. Best collapses are 
achieved with v±_ = 0.49(4). 
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